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Speckle structure of parametric down conversion light has recently received a large attention 
due to relevance in view of applications to quantum imaging 

The possibility of tailoring the speckle size by acting on the pump properties is an interest- 
ing tool for the applications to quantum imaging and in particular to the detection of weak 
object under shot-noise limit. 

Here we present a systematic detailed experimental study of the speckle structure produced 
in type II PDC with particular attention to its variation with pump beam properties. 
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1. Introduction 

Thermal or pseudo-thermal light (as the one obtained by scattering of coherent light 
by a diffuser) presents a random intensity distribution known as speckle pattern 
(J7|). This structure can have interesting applications, e.g. in metrology @. 

In particular, speckle structure of parametric down conversion (PDC) light has 
recently received a large attention due to relevance in view of applications to quan- 
tum imaging @. 

The aim of Sub Shot Noise (SSN) quantum imaging is to obtain the image of 
a weak absorbing object with a level of noise below the minimum threshold that 
is unavoidable in the classical framework of light detection. Being interested in 
measuring an image, one is forced to consider a multi-mode source, which is able 
to display quantum correlation also in the spatial domain. Theoretically, this goal 
can be achieved by exploiting the quantum correlation in the photon number be- 
tween symmetrical modes of SPDC . Typically the far field emission is collected 
by a high quantum efficiency CCD camera. It is fundamental to set the dimension 
of the modes coherence areas with respect to the dimension of the pixels. In par- 
ticular, the single pixel dimension must be of the same order of magnitude of the 
coherence area or bigger, in order to fulfill the sub-shot noise correlation condition, 
compatibly with the requirement of large photon number operation. Thus, the pos- 
sibility of tailoring the speckle size by acting on the intensity and size of the pump 
beam represents an interesting tool for the applications to quantum imaging and 
in particular to the detection of weak objects under shot-noise limit ([3J 0- 
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A detailed theory of correlations and speckle structure in PDC has been devel- 
oped in Q and, in another regime, in (0) furthermore, experimental results were 
presented in 

Nevertheless, a systematic comparison of experimental variation of speckle size 
and their correlations with the theoretical results of @ is still missing. 

In this paper we present a systematic detailed experimental study of the speckle 
structure produced in type II PDC with particular attention to its variation with 
pump beam properties. In particular dependence on pump power and size are 
investigated in detail: results that will represent a test bench for the theoretical 
models. 



2. Theory 

The process of SPDC is particularly suitable for studying the spatial quantum 
correlations because it takes place with a large bandwidth in the spatial 

frequency domain. Any pair of transverse modes of the radiation (usually dubbed 
idler and signal), characterized by two opposite transverse momenta q and — q, are 
correlated in the photon number, i.e. they contain, in an ideal situation, the same 
number of photons. In the far field zone, the single transverse mode is characterized 
by a coherence area, namely the uncertainty on the emission angle $ (tani? = 
Xq/2Tr, A being the wavelength) of the twin photons. It derives from two effects 
that participate in the relaxation of the phase matching condition. On the one 
side the finite transverse dimension of the gain area inside the crystal, coinciding 
with the pump radius w p at low parametric gain. On the other side the finite 
longitudinal dimension of the system, i.e. along the pump propagation direction, 
that is generally given by the crystal length I. 

The appearance of the emission is a speckled structure in which the speckles 
have, roughly, the dimension of the coherence area and for any speckle at position 
q there exists a symmetrical one in — q with equal intensity. This is rather evident 
in the ccd images of SPDC shown in Fig. [TJ 

In the following we summarize, very briefly, some elements of the theory describ- 
ing this structure. 

Omitting some unessential constants, the Hamiltonian describing the three fields 
parametric interaction is 

Hj(t) oc ^ X (2) Ei +) (r,t) E^\r,t) E^(v,t) d 3 r + h.c (1) 

The evolution of the quantum system guided by Hamiltonian ([I]) , in the case of 
relatively high gain regime and non-plane- wave pump, requires a numerical solution 
and it is discussed in detail in Q. Anyway, in the low gain regime an analytical 
solution is achievable Therefore, it is worth to briefly mention the result 

in the first order of the perturbation theory {g <C 1) for a gaussian pump, where 
the quantum state of the scattered light has the entangled form 



1 see also Q for the seeded case. 
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(a) 




(b) 




Figure 1. CCD's images of the speckled structure of the type II PDC emission from a BBO non-linear 
crystal. The two partial rings in (a) represent the two correlated beams 1 and 2, selected by a narrow filter 
(Bandwidth=10 nm) around the degeneracy (uji = u>2 = lu p /2), whereas in (b) no filtering is applied. On 
the right hand site, are shown the respective cross-correlation functions, evaluated by fixing a region Ri in 
the signal field and moving a region R2 in the idler field. A point in the cross-correlation graph correspond 
to the correlation coefficient between Ri and R2 for a single position of i?2- The axes report the pixel by 
pixel displacement of region R2- 



I vac) + exp 



Hidt 



|0> 



|vac) + ^2 "52 F (<lii<l2,ty liqi-n) 1^2 > 
qi,q2 U 



(2) 



F(q 1 ,q 2 ,n) = g ■ sine 



Afc(qi,q 2) Q) • I 



-(qi+q 2 ) 2 



uj\ = ujp/2 + ri, ui2 = ojp/2 — n. 



(3) 



qi 5 2 being the transverse wave vectors and u p ,ui,u}2 the pump, idler and sig- 
nal frequencies respectively. The coherence area, in the limit of low parametric 
gain g, can be estimated by the angular structure of the coincidence probability 
|i ? (qi,q2)| 2 at some fixed frequency 0. As mentioned before, now is clear that 
we deal with two functions that enter in the shaping of the coherence area: the 
sine function and the Fourier transformed gaussian pump profile. Since they are 
multiplied, the narrower determines the dimension of the area. The Half Width 
Half Maximum of the gaussian function, appearing in (ph, is 5q = y/2\n(2)/w p . 
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If we expand the longitudinal wave detuning around the exact matching point 
Afc(qo,qo,^), the linear part (J77J) dominates for angles $o(qo)) n °t too next to 
collinear regime and the sine function turns out to have a HWHM of Aq = 

2, 78/(1 tan $o) at degeneracy. On the other hand, around the collinear emission, 
the quadratic term prevails @ and the bandwidth becomes Aq = 2, 78*(27r/AZ) 1 / 2 . 
Concerning our experiment, we consider small (but not zero) emission angles $ and 
large enough pump radius w p , such that we always work in the region Sq/Aq < 1. 
Therefore, in principle, the dimension of the coherence area is only determined by 
the pump waist. 

When moving to higher gain regime, which is of greater interest for our experi- 
ment, the number of photon pairs generated in the single mode increases exponen- 
tially as oc sinh 2 (<7) i.e. a large number of photons is emitted in the coherence time 
along the direction In this case, also the pump amplitude becomes important 
in the determination of the speckles dimension. As described in @, this can be 
explained by a qualitative argumentation: inside the crystal, the cascading effect 
that causes the exponential growth of the number of generated photons is enhanced 
in the region where the pump field takes its highest value, i.e. close to the center 
of the beam. Thus, in high gain regime, most of the photon pairs are produced 
where the pump field is closed to its peak value. As a result the effective region 
of amplification inside the crystal becomes narrower than the beam profile. Thus, 
in the far field one should consider the speckles as the Fourier transform of the 
effective gain profile, that being narrower, produces larger speckles. 

A further fundamental consideration for the practical implementation is that in 
high gain regime, instead of measuring the coincidences between two photons by 
means of two single photo-detectors, one collects a large portion of the emission 
by using for instance a CCD array with a certain fixed exposure time. Within this 
time several photons are collected by the single pixels and the result is an intensity 
pattern, having the spatial resolution of the pixel. Looking at the images, we can 
appreciate the speckled structure and a certain level of correlation of the speckles 
intensity between the signal and idler arms (Fig. ga) and (b)). We can define 
also the auto-correlation function of the signal intensity pattern itself, since in the 
single transverse mode of the signal arm there are many photons. To be precise the 
speckle's dimension is better related to the spread of this function, although the 
two functions, the cross- and the auto- correlation, present in the very high gain 
regime the same behaviour with respect to the pump parameters, see Q. 

From the experimental view-point it is convenient to study the auto-correlation 
because of the higher visibility that allows a more accurate estimation of its size. 

3. Experiment 

Our setup is depicted in Fig. |3j A type II BBO non-linear crystal (1 = 1 cm) is 
pumped by the third harmonic (wavelength of 355 nm ) of a Q-switched Nd:Yag 
laser. The pulses have a duration of 5ns with a repetition rate of 10 Hz and a 
maximum energy, at the selected wavelength, of about 200 mJ. The pump beam 
crosses a spatial filter (a lens with /=50 cm and an iris of 250 /mi of diameter), 
in order to eliminate the non-gaussian components (see fig. [2]). Before entering 
the crystal the pump is collimated and its diameter is varied, when necessary, by 
changing the distance between two lenses (a biconvex and a biconcave) placed 
after the spatial filter. After the crystal, the pump is stopped by a UV mirror, 
transparent to the visible (T=87% measured at 685nm), and by a low frequency- 
pass filter. The down converted photons (signal and idler) pass through a lens of 5 
cm of diameter (/ = 10 cm) and an interference filter centered at the degeneracy 
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Figure 2. Image of the pump beam before (on the left) and after (on the right) spatial filtering. 

A=710 nm (lOnm bandwidth) and finally measured by a CCD camera. We used 
a 1340X400 CCD array, Princeton Pixis:400BR (pixel size of 20 /mi), with high 
quantum efficiency (80%) and low noise (A = 4 electrons/pixel). The far field 
is observed at the focal plane of the lens in a / — / optical configuration, that 
ensures that we image the Fourier transform of the crystal exit surface. Therefore 
a single transverse wavevector q is associated to a single point x = (A//27r)q in 
the detection plane. The CCD acquisition time is set to 90 ms, so that each frame 
corresponds to the PDC generated by a single shot of the laser. 

Let us define JVr(x) the intensity level, proportional to the number of pho- 
tons, registered by the pixel in the position x of a region R. 5N R (x) = N R (x.) — 
(N R (x)} is the fluctuation around the mean value that is estimated as (N R (-x)) = 
(1/n) ^x-^k( x )> with n the number of pixels. We evaluate the normalized spatial 
auto-correlations of the intensity fluctuations by choosing a large arbitrary region 
R, belonging for instance to the signal portion of the image, and measuring 

y/{5N R (x)*) (5N R (x + 2 )' 

where £ is the displacement vector in the pixel space and assumes discrete values. 
C(£) has unit value for £ = 0. 

Fig. [4] shows some typical section of the auto-correlation for different value of 
the pump intensity. One note that, especially for small photon number, the auto- 
correlation is dominated by the peak at £ = and drop down even for displacement 
of one pixel, creating symmetrically a sort of shoulder. This is explained by the 
presence of noise at the scale of the pixel: first of all the shot-noise that gives a 
contribution of (N R (x)), then the electronic noise in the read-out and digitalization 
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Figure 3. Experimental setup. A triplicated Nd-Yag laser beam, after spatial filtering, produces type II 
PDC in a BBO crystal, which is then measured, after an interference filter and pump elimination, by a 
CCD camera. 
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Figure 4. Section of the two dimensional auto-correlation function evaluated in single shot images at 
different pump intensity 



process is A and, finally, the quantum efficiency that fluctuates with a standard 
deviation 5rj = 3% (we evaluated it using a self-calibration procedure similar to 
the one described in (jjgp ). On the other side the fluctuations at the scale of the 
coherence area of the light are given by (-/V#(x)) 2 /M, where M is the number 
of modes (spatial modes * temporal modes), collected by one pixel. Therefore a 
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large number of spatial modes limits the visibility of the genuine auto-correlation 
of the field with respect to the other sources of noise due to the detection by CCD. 
Anyway we can take the width of the shoulder as the most correct estimation of 
the coherence radius. Fig. [5] reports on the left hand side, one of our main results, 
namely the speckle's radius R co h as a function of the photon number per pixel for 
two different pump size w p = 1.0 mm and w p = 0.7 mm. The increasing of the 
radius of the auto-correlation function versus the mean number of photons can not 
be explained by the theory in the low gain regime, predicting no dependence on the 
gain, i.e. on the pump intensity. This is a signature of the fact that we are working 
in a non linear regime of PDC in which the parametric gain g does not fulfill the 
condition g <C 1, as we will discuss in the next paragraph. It must be noticed that 
we are constrained in the range of intensities of the pump. For high intensities we 
are limited by the damage threshold of optical components, for low intensities by 
the visibility of the speckle structure because we collect a lot of temporal modes 
in the same frame. 

In order to comment the behavior of the radius it is important to evaluate the 
effective gain region. We adopt two approaches. First we try to evaluate g by the 
statistics of the light. Approximatively, according to Fig. [5] we can consider the 
dimension of the pixel always smaller than the coherence area and thus the number 
of collected spatial modes ~ 1 ( \14\l - We evaluate the number of temporal modes M, 
within this approximation, by using the property of the thermal statistics corrected 
for the presence of experimental imperfections described above: 



(5 2 N) = (N) + -^-^ — h 5 2 7]{N) 2 + A 2 (5) 

The estimation is performed on a statistical ensemble of independent pixel belong- 
ing to a single shot frame, thus is definitely a spatial statistics. 
After that, we can approximatively estimate the parametric gain by the relation 

(N) = 7 1 -i lco u-M -smtfig), (6) 

where sinh 2 (g) is the number of photons in the single space-temporal mode, r} is 
the quantum efficiency (including losses) and rj co u is the collection efficiency. It 
takes into account that the pixel collects only a portion of the spatial mode and 
roughly is given by the ratio between the pixel area and the auto-correlation area 
rjcou ~ Api x /A co h. The result is reported in Fig. |5j on the right hand side. As 
expected, the gain g is larger than 1, confirming that we are working in high gain 
regime. Let also notice that, although each point is an average over tens of images 
they are distributed rather noisy, as well as the radius (on the left insert). 

In order to explain the main features of the data it is fundamental to take into 
account the characteristics of the pump. It is the third harmonic of a multimode 
high power laser and the power is varied by the delay between the Q-switch turn-on 
and the lamp flash. The calibration curve delay-power (mean pulse power) has been 
measured by a power meter, and we observed a reproducibility with uncertainty 
around 10%. The fluctuation of the mean pulse power is about 20% (measured by 
a photodiode). Since the mean number of photons is proportional to sinh 2 (g), the 
fluctuations of the pump power generate large fluctuations in the photons number 
from frame-to-frame. 

Furthermore, a multimode laser when the number of modes n is larger than few 
units should present a thermal statistics of the intensity fluctuation according to 
the formula 51 /I = \J\ — 1/n fli^p and, after the non linear process of second and 



:11 Journal of Modern Optics speckle 



8 







* Coh. Radius [pix] (Rpump=0.50 mm) 

* Coh. Radius [pix](Rpurnp=0.35 mm) 










• 

• 

«• 


• 


• 

1 




♦ 







Gain by Therm. Statistics (Rpump=0 50 mm) 
• Gain by Therm Statistics (Rpump=0 35 mm) 



+ * * 



1000 20D0 3000 4000 50D0 6000 



../ 

* : 



100D 2000 3000 4000 5000 6000 



Figure 5. On the left, speckle's radius and, on the right, parametric gain estimation, both as function of 
photon number per pixel, for two pump radii 




Figure 6. (a) Number of photons in a single mode as function of pump amplitude (b) Parametric gain as 
function of the photon number per pixel estimated from (a) and from the property of the thermal statistics 
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Figure 7. Coherence radius as function of the parametric gain 



third harmonic generation the fluctuations can increase even more. The expected 
temporal profile of the output pulse is very far from a smooth gaussian function. 
Therefore, taken into account the strong non-linearity of the PDC generation, it is 
reasonable to consider that, inside the single pulse entering the crystal, only few 
peaks of intensity contribute to the large majority of the PDC emission. The main 
consequence is that the number of effective temporal modes M, collected in the 
single shot, is much less than what expected by evaluating the ratio between the 
pulse duration (5 ns) and the PDC coherence time r p ~ 1 ps and M also fluctuates 
randomly from image to image. By using Eq. [5] we obtain a mean number of 
effective temporal modes M = 170 ± 80 when w p = 1.0 mm and M = 60 ± 40 
for w p = 0.7 mm. At the same time, the values of the parametric gain g are in 
general bigger than 1 (see Fig [5| . The gain region in which we are working ensures 
the possibility to reach a sufficient non-linear cascading effect in the photon pairs 
production inside the crystal. 

As a confirmation of the high non-linearity obtained in our experiment we esti- 
mate the gain g in an independent way In Fig. |6^a) we present the measure- 
ment of the number of photon in the spatial mode defined as (N mo d e ) = (N)-A co h by 
varying the mean pump amplitude A pump , while the pump size is fixed to w p = 1.0 
mm. Here the mean pump amplitude is the square root of the average power of 
the single pulse divide by the transverse pump area measured by a CCD. The data 
are fitted by the equation (Nr) = k ■ smh. 2 (aA pump ), where k is fixed to k = rjM, 
while a (crA pump = g) is the free parameter. The experimental values, mediated on 
three different acquisitions, are a = 2.53 ± 0.04[mm 2 /W] 1 / 2 . Fig. gb) shows the 
comparison between the values of the gain obtained in this way (green squares) 
and in the previous statistical approach (the red dots). The data are very close, 
and we can conclude that we are actually in a strong non linear regime. 

Now, let us consider the trend of the coherence radius versus the parametric gain 
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Figure 8. Observed dependence of the radius (in pixels) of the speckles from the pump diameter 



for the two different pump size as reported in Fig. [7] The data are well fitted by 
two straight lines (not showed in figure), Ri(g) = aig + Ri(0) with i = 1,2. We 
obtain that ai(0.67± 0.06) and a2(0.55 ± 0.06) agree within the their uncertainty. 
Although a comparison with a theoretical model in high gain regime is needed for 
a proper interpretation of our data, we can consider this behavior as an evidence 
of the influence of the gain and of the pump diameter on the dimension of the 
speckles. In fact, for the same value of the gain, the speckles radius is always larger 
when the pump size is smaller as predicted in the low gain. 

Finally, we investigate the dependence of the radius of the speckles from the pump 
diameter, as shown in figure [8| We fix the power of the laser to 0,78 MW, within 
the fluctuations, and we change the diameter varying the distance between the two 
collimating lenses (see figure [3]). Since the diameter changes, the intensity and the 
gain g change as well. The theory, in low gain regime, provides that the radius of 
the speckles is proportional to the inverse of the pump size (w p ) obtaining the red 
line in figure. It approach the experimental data only when the pump diameter is 
large. This confirms the role of the high gain regime in the speckle size. In fact, 
together with the reduction of the pump diameter, the gain increases, and thus the 
effective gain area is even more reduced. This effect impresses upon the speckles 
size a stronger dependence with respect to the pump size. 



Conclusion 



This paper provides a detailed experimental study of the size of the coherence area 
in PDC in the high gain regime. We show that the speckles present, not only a 
dependence on the pump diameter, as in the usual low gain regime, but also a 
strong dependence form the pump intensity. The understanding of the behavior 
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of the coherence area in the high gain regime is fundamental for the innovative 
application in the field of quantum imaging. Our results provide the basis for a 
comparison with the theoretical models; in general, the observed dependencies 
of speckles size on pump parameters allow the tailoring of speckles properties, 
an instrument useful not only for applications to quantum imaging but also for 
quantum metrology, quantum information, etc. 
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